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ABSTRACT 

Extremely luminous, "super-Chandrasekhar" (SC) Type la Supernovae (SNe la) are as yet 
an unexplained phenomenon. We analyse a well-observed SN of this class, SN 2009dc, by 
_C ■ modelling its photospheric spectra with a spectral synthesis code, using the technique of 

Q^f "Abundance Tomography". We present spectral models based on different density profiles, 

corresponding to different explosion scenarios, and discuss their consistency. First, we use a 
density structure of a simulated explosion of a 2 Mo rotating C-O white dwarf, which is often 
proposed as a possibility to explain SC SNe la. Then, we test a density profile empirically 
inferred from the evolution of line velocities (blueshifts). This model may be interpreted as 
a core-collapse SN with an ejecta mass of ^3 M©. Finally, we calculate spectra assuming 
an "interaction scenario". In such a scenario, SN 2009dc would be a standard WD explosion 
with a normal intrinsic luminosity, and this luminosity would be augmented by interaction of 
the ejecta with a H-/He-poor circumstellar medium. We find that none of the models tested 
easily explains SN 2009dc. With the 2 M Q WD model, our abundance analysis predicts small 
amounts of burning products in the intermediate-/high-velocity part of the ejecta (v > 9000 
km s _1 ). However, in the original explosion simulations, where the nuclear energy release per 
unit mass is large, burned material is present at high velocities. This contradiction can only 
be resolved if asymmetries strongly affect the radiative transfer or if C-O white dwarfs with 
masses significantly above 2 M Q exist. In a core-collapse scenario, low velocities of Fe-group 
elements are expected, but the abundance stratification in SN 2009dc seems "SN la-like". The 
interaction-based model looks promising, and we have some speculations on possible progen- 
itor configurations. However, radiation-hydrodynamics simulations will be needed to judge 
whether this scenario is realistic at all. 

Key words: supernovae: individual (SN 2009dc) - supernovae: general - techniques: spec- 
troscopic - radiative transfer 



1 INTRODUCTION 



Type la supernovae (SNe la) are subject to vigorous research ow- 
ing to their relevance to different fields of astrophysics. Their pro- 
duction via different paths of binary evolution is investigated in 
detail (e.g. iRuiter et alj|2009l : Idaevs et alj|20ld: IPodsiadlowskil 
20ld). as is their impact on galaxies (e.g. IScannapieco et alj 2008; 



Tang etalj|200^) . Most prominently they are u sed today in cos 



mological studies (e.g . 



Perlmutter et al 



IPerlmutter etaill 19991 : lAstier et alJbOOoHWood-Vasev et alif^ 



1997; Riess et al. 199: 



to measure luminosity distances. SNe la show differences in their 
intrinsic luminosity. Observations have, however, shown that the lu- 
min osity of each object can be estimated from its light curve shape 
(e.g. |Pmllipd[l99ll ; |Phillips et al.lll999h . which is measurable in- 



dependently of distance. Empirical calibration procedures exploit 
this, allowing for the construction of accurate SN Hubble diagrams. 

The fact that SNe la constitute a relatively homogeneous 
SN subclass is usually explained by the homogeneity of their 
progenitor systems, which are supposedly carbon-oxygen white 
dwarfs (WDs) in binary systems. In the "classical" single- 
degenerate Chandrasekhar-mass scenario, one WD explodes af- 
ter having accreted matter from a companion star. A thermonu- 
clear runaway occurs when the W D is close to the stability 
limit faillebrandt & NiemeveJboOOh . i.e. its Chandrasekhar mass 
(Mch.non-iot ~ 1.4 Mq without rotation). Numerical simulations of 
such explosions have been successful in explaining the "leading- 
order" homogeneity in th e observations, and also part of the di - 
versity among the objects jMazzali et al]|2007l ; iKasen et alj|2009h . 
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The most important variation in the properties of SNe la is caused 
by a different efficiency in producing 56 Ni, which powers the op- 
tical light curve through the radioactive decays Ni — > 56 Co — > 
56 Fe. 

Despite the success of standard explosion models, peculiar 
SNe l a, discovered in in c reasing numbers (e.g. iFilippenko et all 
1 1991 IPhiliipsetalJl 19921 : Li et aljEjol IHowell et alj|2006h . are 
difficult to explain thus far. Meanwhile, the precision of cosmo- 
logical studies based on SNe la has significantly increased. In 
this situation, we have to clarify what mechanisms can produce 
SNe la in order to avoid errors in SN cosmology which could 
be ca used by populations of peculiar SNe (cf. e.g. iFolev et all 
l20ld) . Within the "classical" single-degenerate Chandrasekhar- 
mass scenario, there is the possibility of having an "supermas- 
sive" white dwarf as p r ogenitor, which is stabi l ised by rotation (e.g . 
Steinmetz etafl 1 19921: lYoon & Langeil 120051: IHowell et alj 12006 



Pfannes et al J 1 20 1 Oal ; iFinkl 1 20 1 (j : lHachisu et alj|20ll Fink et al. 
2012 in prep.). Also other explosion scenarios are re-considered 
to produce peculiar and even normal SNe l a - most prominently , 
the double-degenerate/WD-merger scenario l lPakmor et alj 201 Oal) 
and the single-degenerate su b-Chandrasekhar scenario i Fink et al] 
l2007ll20ld : lsim et alj2O10h , where a He shell around the WD det- 
onates and triggers the explosion of the star. Furthermo re, there are 
relati vely new, somewhat "exotic" SN la models (e.g. iBravo etlrtl 
l2009h . 

On the observational side, one prominent SN la subclass 
waiting to be explain ed are the "super-Ch a ndrasekhar m a ss" ob - 
jects [SC SNe la - IHowell et alj d2006l) : Ifficken et alj d2007l) : 
IScalzo et al. and references below; the designation refers 

to the Chandrasekhar mass for non-rotating WDs]. The early-time 
spectra of these SNe resemble thos e of other SNe la w hile their lu- 
minosity is abou t twice the ave rage dHowell et al .120061) . According 
to Arnett's rule jArnettJll982h . Ni masses of more than a Chan- 
drasekhar mass are required to explain the luminosity of some of 
these objects. A lu minosity prediction f rom the light-curve shape, 
using the relation o f|PhillipsetalJdl999l) . fails at least for some SC 
SNe dTaubenberger et alj|201 lh . Possibly, part of the luminosity of 
SC SNe comes from energy sources other than 56 Ni decay, such 
as interaction of the ejecta with circumstellar material. The more 
conventional paradigm, however, is that SC SNe la are more mas- 
sive than normal objects (possibly the result of a rotating progenitor 
WD), and pro duce more 56 Ni tha n a Chandrasekhar-mass explo- 
sion could do l lHowell et al1l2006l) . In any case, in-depth analyses 
of observations with radiative transfer models are needed to better 
understand SC objects. 

In this work, we study the recent SC SN 2009dc, which 
reached a peak Mt> i ~ —20 and exploded in a tidal bridge be- 
tween two galaxies (UGC 10063/64, types SO/SBd). Observations 
with an extraordinarily good time coverage are available for this SN 
dYamanaka et alj|2009l: hanaka et alj [20101 : ISilverman etal]|201ll ; 
iTaubenberger et alj201 lh . At maxi mum light, SN 2009d c has spec- 
tral properties similar to SN 2003fg dHowell et al. 1 20061) . and there- 
fore is representative at least for part of the SC objects. Other 
SC SNe are differ ent to some degree. For example, SN 2007if 
dScalzoetal]|2010h shows higher-temperatur e spectra and some - 
what higher line velocities (as also SN 2006gz. lHicken et alj2007l) . 
The reasons for this diversity are unclear, as is the theoretical in- 
terpretation of SC SNe. We use a radiative-transfer code to model 
the observed spectra of SN 2009dc, inferring ejecta densities and 
abundances. The method we employ (Abundance Tomography, 
IStehle et al.ll2005l) has pr oved suitable for inferring the properties 
of quite diverse SNe (e.g. lTanaka et alJl201ll.lMazzali et alj|2008l. 



lHachinger etai]|2009l) from data taken in the photospheric phase, 
when the core of the ejecta is opaque to optical radiation. Assuming 
different density profiles and rise times, we infer a best-matching 
abundance structure, respectively, by fitting the observations with 
synthetic spectra. Comparing models based on different assump- 
tions, we assess the applicability of different explosion scenarios 
by judging the fit quality and the physical consistency of the mod- 
els. 

Below (Sec.O, we first give a relatively extended overview 
of our methodology as it is relevant to the present work. We then 
show a tomography model based on a density profile representing 
the explosion of a 2.0 Mq rotating white dwarf, and another model 
based on an empirically constructed 3.0 Mq density profile (Sec- 
tions [3] and @. Finally, we present a third tomography experiment 
with a more speculative set-up: we model the spectra assuming 
that the SN is a Chandrasekhar-mass explosion (of a non-rotating 
WD) where some light is "externally" generated through interac- 
tion of the ejecta with a surrounding medium (Sec.O- We discuss 
all our models in terms of consistency and theoretical interpretation 
(Sec.[6]( before we conclude (Sec. [7}. 



2 METHOD AND DATA 
2.1 Radiative-transfer code 

We use a 1-D stationary Monte-Carlo radiative-transfer code 



dAbbott & Luc\Jll985|.lMazzali & Lucvlll993l ILucvII 19991. iMazzalil 
l200d and lStehle et al.ll2005h to calculate synthetic SN spectra from 
a given density and abundance profile. 

The code computes the radiative transfer from an assumed 
photosphere through the SN ejecta at a given epoch. An initial den- 
sity profile is used to calculate the densities within the ejecta, as- 
suming homologous, force-free expansion. This expansio n phase 
comm ences shortly after explosion in most SN models le.g. lRopkel 
d2005l) : cf. however Sec. 15. II . so that r = v x t holds for each 
particle at all times relevant to this work in good approximation (r: 
distance from the centre, t: time from explosion, v: velocity at the 
onset of homologous expansion). Radius and velocity can then be 
used interchangeably as spatial coordinates. 

From the photosphere, which is located at an adjustable v p h, 
continuous black body radiation [TJ~ = B v (T p h), where I„ is the 
specific intensity in outward directions, and B v (T v h) the specific 
intensity according to the Planck law at temperature T P h] is as- 
sumed to be emitted into the atmosphere. Notable deviations from 
this approximation to the flux in the inner layers appear in the red 
and infrared, where a discrepancy in the flux level of synthetic to 
observed spectra therefore is sometimes unavoidable. 

The photon packets simulated undergo Thomson scattering 
and interactions with lines, which are treated in the Sobolev ap- 
proximation. A downward branching scheme ensures a good ap- 
proximation to the bound-b ound emissi vity. Radiative equilibrium 
is enforced by construction dLucv|[l999T) of the Monte-Carlo simu- 
lation. The excitation and ionisation state of the matter is calculated 
from the radiation field statistics using a modi fied nebular approx- 
imation dMazzali & Lucvlll993l : lMazzalill2000l) . In this approxima- 
tion, the gas state in each radial grid cell is mostly determined from 
a radiation temperature Tr and a dilution factor W. Tr corresponds 
to the mean frequency of the radiation field and W parametrises its 
energy density (for given Tr). We iterate excitation, ionisation and 
the radiation field in turn until the Tr values within the atmosphere 
are converged to the per cent level. Within these iterations, T P h is 
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automatically adjusted in order to match a given bolometric out- 
put luminosity Lboi. Thus, the photospheric luminosity is adapted 
so as to compensate for reabsorption of radiation, which occurs 
when packets re-enter the photosphere. After the temperature itera- 
tions, the emergent s pectrum is c alculated by formal integration of 
the transfer equation ( lLucv|[l999h , using a source function derived 
from the Monte Carlo statistics. 

2.2 Spectral models 

To summarize, a "SN model" from which our code calculates a 
synthetic spectrum is specified by the density profile, the abundance 
stratification, the luminosity Lboi, the photospheric velocity v P h and 
the time from explosion t. 

Spectral modelling means adjusting this model to obtain the 
best possible match of synthetic and observed spectrum; the model 
parameters will then represent the SN properties. Usually, a spe- 
cific density profile and rise time is assumed before starting to fit 
the spectra. Below, we briefly discuss how the various remaining 
parameters influence the radiative transfer and the spectra, and we 
describe our general strategy for fitting the observations, including 
the tomography technique. 

2.2.1 Abundances and their imprint on the spectrum 

The part of the abundance structure relevant for the SN spectrum 
varies with the epoch (see also Sec. l2.2.5l >. For each model, only the 
abundances above the photosphere (i.e. at v > v p h) are relevant. 

The abundances of intermediate-mass elements (IME, e.g. 
Mg, Si, S, Ca) normally influence the strength of individual spectral 
features. In contrast, most Fe-group elements do not leave distinct 
features in the optical spectra. They rather influence the spectrum 
by modulating the line blocking in the UV and by a "transfer" of 
flux to the red [via branching, see iMazzalil J200(ih l. The line block- 
ing effect has the consequence of backwarming within the atmo- 
sphere (steepening of the temperature profile), exerting an indirect 
influence on line strengths of all elements via the ionisation and 
excitation state of the plasma. This state is not fixed by input pa- 
rameters in our approach, but is determined such that it is consistent 
with the radiation field. 

2.2.2 Influence of Lboi and v p h 

Lboi is the parameter controlling the flux level in the synthetic spec- 
trum. For each epoch, it is adapted during the fitting process so as 
to match the overall flux in the optical part of the observed (and 
flux-calibrated) spectrum. Therefore, it is easily constrained to a 
reasonable precision. 

The photospheric velocity v p h has various, more complicated 
effects on the model. First, it sets a lower limit to the radii at which 
line absorption can occur in the simulation. Second, it determines 
the integrated optical depth of the atmosphere, which should be ap- 
propriate (r > | in the relevant wavelength range) even if the no- 
tion of a sharp photos phere emitting a bl ack body is quite approxi- 
mative in Type I SNe (Sauer et al. 200^). For our simulations, it is 
important that the photosphere is assumed deep enough such that 
the assumption of a black-body like emission at that depth is satis- 
fied as well as possible. A very low v p h, however, generally leads 
to a very blue radiation field at the lower boundary, as the emitting 
surface shrinks. When this is not fully compensated by attenuation 
above the photosphere, the model spectrum will be bluer. Within 



the freedom these constraints give us (roughly ± 1000km s _1 ), we 
choose Vph such that we obtain an optimum fit. 

During the fitting process, we assess and take into account 
additional effects of Lboi and v p h on the strength of spectral lines 
via the plasma state. A lower v p h tends to make the radiation field 
bluer at least in parts of the atmosphere (cf. above). This results 
in higher radiation temperatures and an increased ionisation (and 
occupation of excited levels). A higher luminosity Lboi makes the 
radiation field stronger and bluer, with the same result. 



2.2.3 Time t - sensitivity of the model on the rise time 

The time from explosion t must be determined from the respective 
observational epoch [usually given relative to B-band maximum 
(B max.)] and the assumeqj rise time t T , in which the SN rises to 
B max. after explosion. 

When the density distribution in velocity space is kept con- 
stant, and t t is changed, the spectrum shows an overall shift of 
line velocities (mean blueshifts) for the following reasons: 

• Models with a smaller t (smaller t r ) usually have photospheres 
at higher v v h. With the more compact ejecta, the higher v P h is 
necessary in order to avoid excessive backscattering and a too 
blue spectrum of the outgoing radiation at the lower boundary 
(when a high luminosity is emitted from a small photosphere). 
The higher v p h then leads to some increase of line velocities. 

• The ejecta are more compact for a smaller t r , since 
r — v x t , where t is the sum of t r and the epoch (relative 
to B maximum). Therefore, for a smaller t r the energy den- 
sity of the radiation field inside the ejecta tends to be larger at 
any given velocity v (with Lboi kept constant). The zone where 
ionisation is low enough for line formation (i.e. triply- and es- 
pecially doubly-ionised species dominate) will be located at 
higher v. This then reflects in higher line velocities. 

When reducing t, below realistic values, line formation would 
have to occur in the outermost zones of the ejecta. The low densities 
(and the strong ionisation due to the low densities) in these zones 
however disfavour the formation of prominent features such as Si II 
A6355. As a result, the spectra of an observed object can not be 
explained if too low a rise time is assumed. 



2.2.4 Fitting observed spectra 

The process of spectral fitting begins with choosing a density pro- 
file and a rise time, which subsequently remain unchanged. We then 
iteratively determine optimum values of those abundances which 
influence features on the one hand, and parameters that control the 
temperature structure and the line velocities on the other. We try to 
match the spectral features as well as the overall flux, but if this is 
not possible simultaneously, we prefer a better fit to the lines. In 
the present case, we sometimes allow the synthetic spectra to be 
too red if otherwise the model atmospheres would become too hot. 



1 The rise time t r can onl y to some degree be constra ined from the available 
early light curve data [cf.Taubenbeiger et al. ( 20 11) and Sec. |2.4| . 

2 Here, we mean that the time-invariant density distribution ^j^(v), i.e. 
the explosion model, is kept constant. 
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2.2.5 Abundance Tomography 

The method of Abundance Tomography iStehle et al. 120051) ex- 
tends the approach of fitting one spectrum dMazzali&Lucvlll993h . 
It aims at inferring the abundance profile of a SN step by step from 
the outer to the inner layers (assuming a density structure and a rise 
time). To this end, one fits a time series of photospheric spectra, in 
which different layers of the ejecta leave their imprint. In the pho- 
tospheric phase, some core of the ejecta is optically thick at UV 
and optical wavelengths. A spectrum is then only sensitive to the 
optically thin part of the ejecta. As time progresses, a larger and 
larger fraction of the ejecta becomes visible. Therefore, later spec- 
tra are sensitive to more central layers, while the influence of the 
outer layers on them is smaller. 

In order to create a model for the outermost ejecta of an ob- 
ject, we fit the earliest spectrum available. The photospheric veloc- 
ity, the luminosity and the abundances are optimised. In the case of 
SN 2009dc, we use three abundance zones above the earliest pho- 
tosphere in order to have composition gradients in the outer ejecta, 
allowing for an optimum fit (cf. Sec. l3.lt . 

The spectrum at the next epoch carries the imprint of the mate- 
rial in the outer envelope and additionally that of the layers which 
the photosphere has receded into in the meantime. As the outer 
abundance structure is already constrained, the new photospheric 
velocity and the abundances of the newly-visible layers are now 
the primary fitting parameters. If the abundances of the outer layers 
have a significant impact on the quality of the fit, and no decent 
fit is possible with the values inferred earlier, those values are re- 
vised (normally by less than 30% of their initial value) to equally 
optimise the fit to all observations. The procedure is continued with 
later spectra until the complete time series is fitted. 

Degeneracies, i.e. equally good fits with different parameters, 
are the most important source of uncertainty in the fitting process. 
Abundances not well constrained by an early spectrum can some- 
times be determined when fitting a later spectrum. The respec- 
tive values may then undergo significant revision (> 30%). In a 
few cases, this a-posteriori determination does not work. Then, we 
assume the unknown abundances to match typical nucleosynthe - 
sis patterns within SN la explosion models l llwamoto et al.|[l999h . 
Most importantly, when Fe-group elements do not lead to the for- 
mation of individual features, the mix of Fe-group elements is set 
so as to resemble that in the outer, intermediate or inner zone of the 
llwamoto et al. I dl999l) models as appropriate. Furthermore, in the 
intermediate zone of the ejecta, Si has been assumed to be the most 
abundant element, even when these zones contributed little to the 
formation of the spectral features. 



2.3 Study of SN 2009dc: concept 

SN 2009dc is in many respects different to other SNe which Abun- 
dance Tomography has been applied to. In particular, as it is un- 
known which explosion model applies, the density profile and the 
rise time t r constitute major elements of uncertainty in our analysis. 

We reflect this in our modelling approach: we perform Abun- 
dance Tomography three times, using different density profiles 
which represent different explosion types. The aim is to compare 
the different models in the end and to assess their consistency. Be- 
fore conducting our final tomography, respectively, we decide on 
an optimum rise time for each density profile. 



2.4 Observations of SN 2009dc 

We analyse the luminosi ty-calibrated spectra of SN 2009dc of 
iTaubenberger et al] d201 ll) . This data set samples SN 2009dc ex- 
tremely well from 9.4 d before B maximum to late times. We do 
not model spectra later than 36.4 d past B maximum. At this epoch 
- late for a SN la - the photosphere recedes deep into the 56 Ni-rich 
zone and our assumption of radiative equilibrium in the simulated 
atmosphere eventually breaks down. 

As the spectra only show a slow evolution of line velocities, 
we expect a very moderate photospheric recession for all our mod- 
els (which is confirmed by our actual calculations). This has an im- 
pact on the tomography method: if the abundance zone between 
two subsequent photospheres is too thin, the zone has little in- 
fluence on the spectrum, which implies a large uncertainty in the 
abundance determination. In order to avoid this, we only model the 
spectra at -9.4 d, -3.7 d, +4.5 d, +12.5 d, +22.6 d and +36.4 d 
with respect to B maximum, leaving out t he other ones in between. 

Following lTaubenberger et all d20 1 ll) . we generally assume an 
extinction towards SN 2009dc of 

E(B - V) tot ~ E(B - V) host + E(B - V)oa ~ 0.17 mag. 

Because of the relatively short distance to the SN, we can add up 
the reddening values without taking redshift into account here. The 
Galactic extinction ta ken int o account is E(B — V)o a i = 0.07 mag 
dSchlegel et al.ll998]) , while lBurstein & Heilesl Jl982h give a lower 
value E(B — l / )oai = 0.04 mag. The reddening due to the hos t 
galaxy carries an error of 0.08 mag iTaubenberger et al] l201ll) . 
Therefore, the lower-limit reddening within the error bars is 

E(B - V)m > 0.06 mag. 

While we generally use the standard value of 0.17 mag, we repeat 
some analyses for a value of 0.06 mag to test the influence on the 
results. 

The rise time of SN 2009dc has been observationally es- 
timated to_be_shorte^Jhan ~30 d and longer than ~22d [we 
follow ITaubenberger et al] d20 1 lh here who give a large r, more 
conservatively-estimated range than I Silverman et au d201 lh l ■ 



3 TOMOGRAPHY ASSUMING AN EXPLOSION OF AN 
SUPERMASSIVE WD 

Our first set of models is based on a density profile representative 
of an explosion of a rotating WD: we adopt the profile resulting 
fro m the "AWD 3det" simulation of Fink et al. [2012, in prep.; see 
al solFinkl(l2010l)l , who re -calculated the 'AWD3 detonation" model 
of IPfannes et al. (2010b) with improved detonation physics. The 
AWD3det model can be regarded prototypical for explosions of 
rotating, massive (2 Mq) white dwarfs. The thermonuclear flame 
in AWD3det propagates as a pure detonation (shock-driven, su- 
personic combustion wave) from the very beginning. A pure de- 
flagration (subsonic combustion flame) in a massive, rotating WD 
would be less likely to explain SN 2009dc, as it does not produce 
much 56 Ni (~0.7 Mq ), and leads to a strongly mixed composition 
jPfannes et al.ll20l"0a) . Delayed-detonation models for these pro- 
genitors, in which the fl ame initially proc eeds as a deflagration and 
later as a detonation (cf. lKhokhlovlll99ll) . produce roughly similar 
outcomes as AWD3det, as the initial defla gration in a rotating WD 
tends to be weak and burns little material dFinkl2010r) . 

In order to conduct a tomography analysis in spherical sym- 
metry as usual, we had to reduce the AWD3det density field to an 
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Figure 1. Density profiles used for Abundance Tomography (Sections [5] 
[4] and |5J, and W7 profile as a reference. The AWD3det model has been 
converted to ID (cf. Sec. [3); the 09dc-exp model is shown with its slope 
adapted to a rise time of 29 d (cf. Sec. 14. Q . 

averaged 1-D density profile by adding up the cell masses into ra- 
dial velocity bins (Fig.QJ. Asp hericities ar e to some degree inherent 
in a rotating explosion model dFinJdl2010h . We keep this in mind as 
a caveat for the discussion of our results. 

3.1 Rise time estimate 

We determine an optimum rise time for our AWD3det-based mod- 
els from the —9.4 d spectrum of SN 2009dc. The earliest observa- 
tion is most useful for this purpose, as it is most sensitive on the 
assumed rise time (the relative difference in t, the time from explo- 
sion onset, is largest at early epochs). 

The abundance structure of the early-time models comprises 
three abundance shells. The set up allows optimisation, but avoids 
an overparametrisation and remains sensible in terms of the chem- 
ical structure: the models have a zone of thickness 600 km s - 
above the photosphere with some Fe-group elements (up to 15% 
by mass; higher fractions are probably not realistic) and high IME 
mass fraction. Above this zone, there is a layer depleted in Fe-group 
elements but still rich in IME, with a thickness of ~2000 km s _1 . It 
causes relatively little backscattering, such that temperature drops 
quickly with radius and becomes favourable for the formation of 
prominent lines due to singly-ionised IME. Still above, the atmo- 
sphere is assumed to consist of a C-O mix containing only traces 
of other elements (mass fractions generally < 10% of those in the 
zone below). The photospheric velocities are selected such that a 
reasonable temperature structure (i.e. a large enough strength of 
Si II lines) results and less than two thirds of the radiation pack- 
ets are re-absorbed due to backscattering. Whenever in doubt, we 
choose the lowest possible photospheric velocity in order not to 
artificially suppress the formation of lines deep inside the ejecta. 

As described in Sec 12.2.31 the line velocities in the synthetic 
spectra change with t r . We set up models for different t r and com- 
pare the positions of th e synthetic and the ob served Si II A6355 
featur^fl (measured as in lHachinger et alfeoOol) . Precisely, the rise 
times we have tested are t t = 21 d, 22.5 d, 25 d, 27.5 d, 30 d, 35 d, 
and 40 d (only the most promising models have been finalised). For 
the four models which show the least mismatch AA in the position 
of the observed and synthetic Si II A6355 feature, t r and AA values 

3 Si II A6355 veloci ties in normal SN e are known to be well fitted by our 
spectral models (e.g.[Stehle et al. 20(3). 



Table 1. Rise time estimates from our spectral fits using the AWD3det 
density profile described in the text, for different values of the reddening 
E(B—V) and metallicity. The second and third columns give, for each den- 
sity profile, the four values of the rise time t r for which optimised models 
have been calculated, and the corresponding i> pn in these models. The fourth 
column gives the wavelength deviation between the synthetic Si II A6355 
feature and the observed one AA = A syn — A„t> s (with A (, s = 6165.6 A). 
From a linear fit to the AA — t, relation, we calculate an optimum rise time 
tr.opt for each density profile. 



p profile t r 
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[d] 
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standard: E(B - V) = 0.17 ma 


g, prog, metallicity Z ~ 


AWD3dct 27.5 


9440 
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30.0 
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-3.7 




35.0 


7120 


20.8 




40.0 


5950 


42.1 


31.5 


E(B -V) = 


0.06 mag, prog, metallicity Z ~ 





AWD3dct 22.5 


10610 


-42.9 




25.0 


9100 


-26.8 




27.5 


7860 


1.9 




30.0 


6850 


17.0 


27.8 


E(B - V) = 


- 0.17 mag, solar prog, metallicity 


AWD3dct 27.5 


10070 


-28.5 




30.0 


9200 


-8.0 




35.0 


7680 


19.1 




40.0 


6560 


38.7 


32.1 



are given in Table Q] These values sample the relation between t Y 
and AA, to which we calculate a linear fit (least-squares method). 
We demand AA = and obtain from the linear-fit relation an opti- 
mum value fr.opt (see again TableQ}. 

Two example models, where one can clearly see the shift in 
the A6355 position with the assumed rise time, are shown in Fig. 
[2] together with the observed spectrum of SN 2009dc. A detailed 
discussion of the SN spectrum and our final optimum model based 
on AWD3det is given in the next Section (Sec. |3.2t . 

We first performed our rise time experiment assuming a to- 
tal extinction of E(B — V) = 0.17 mag and negligible progenitor 
metallicity. The latter assumption, which proved optimum for the 
fit quality, has two main implications when modelling the interme- 
diate and outer ejecta: First, there are no lower limits to the metal 
abundances in the progenitor WD (which we would have to respect 
when determining the optimum abundance structure). Second, no 
stable Fe is produced in the outer layers which undergo incomplete 
Si-burning in the explosion, as a pure C/O mix provides no ex- 
cess neutrons dlwamoto et al .1 1999b . Hence, there is no zone where 
stable Fe is the dominant iron-group element directly after the ex- 
plosion. When the Fe-group mix cannot be determined from the 
spectra, we therefore always assume a mix of Fe-group elements 
which is dominated 56 Ni and its decay products and/or by 52 Cr (cf. 
Sec. l2~2~5l 

Afterwards, we repeated the analysis assuming solar metallic- 
ity and a lower-than-standard reddening E(B — V) to t ~ 0.06 mag 
(see Sec. |2,4t , respectively. This enables us to estimate the influ- 
ence of these assumptions on our study. The respective results are 
shown in Table[T]as well. 

For AWD3det, with standard assumptions, the optimum rise 
time turns out to be 31.5 d (Table[T]l. Evidently, optimum rise times 
from measured AA values will always carry some uncertainty, from 
the measurements as well as the modelling methods. With one step 
in our t Y grid, the shift in line position (Table QJ is, however, rel- 
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3500 4000 4500 5000 5500 6000 6500 7000 



X [A] 

Figure 2. AWD3det-based models for the —9.4 d spectrum of SN 2009dc for negligible progenitor metallicity and E(B — V) = 0.17 mag. We show the 
models for the most extreme t r values listed in Table[T] An optimum model would fit the velocity of the Si II A6355 line (marks; inset). With a very short rise 
time, the velocity (blueshift) of the synthetic feature is too high (blue graph); with a very long rise time, the velocity is too low (orange graph). 



atively large, so that we can crudely estimate the error in our t T 
determination to be equal or smaller than an average step in our t r 
grid (3 d). We thus give the following limit on the rise time: 

ir, 09dc/AWD3cfet > 28.5 d. 

A very similar limit is found for solar progenitor metallic- 
ity. If one assumes a weaker reddening of E(B — V) ~ 0.06 mag, 
the optimum rise times for our models are lower. This is due to 
the lower luminosity, reducing the ionisation and allowing lines to 
form further inside in the ejecta (i.e. line velocities as low as in 
the observations of SN 2009dc can be explained by models with 
relatively short rise times, cf. Sec. 12.2.3b . From Table [7] we con- 
clude that for E(B — V) ~ 0.06 mag our rise time limit would 
be £ r ,09dc/AWD3det > 24.8 d. However, at present we see no com- 
pe lling reason for preferring the low reddening over the one given 
bv lTaubenberger et al .1 d201ll) . 

3.2 Tomography 

The synthetic spectra based on AWD3det are shown in Fig. [3] The 
rise time we have used for the calculations is 30 d. The optimum 
rise time inferred above is a bit l onger, but 30 d is the stri ct up- 
per limit allowed by observations dTaubenberger et al]|201 ll) . The 
code-input parameters of our tomography models are compiled in 
Appendix |B1 

3.2.1 — 9.4 d and —3.7 d models 

The observed —9.4 d spectrum of SN 2009dc extends down to 
3000 A, where the flux is unusually high for a SN la. Despite the 
high luminosity and blue colour, the optical spectrum shows strong 
signatures of singly-ionised Fe-group and IME species. The simi- 
larly lumin ous SN 2007if does not show clear Si II or S II lines at 
this epoch (Scal zo et al.ll2O10h . This points towards smaller IME 
abundances or perhaps tow ards a higher degree of ionisation in 
S N 2007if [cf SN 1 99 IT - iMazzali et all l l 19951) vs. SN 1999ee 
- lHamuv et al] J2002h , where variations in ionisation are crucial]. 

In the case of SN 2009dc, an optimum fit to the spectra is 
achieved with three abundance zones in the outer ejecta (just as 
described in Section [3~"Q . From the outside inwards, these zones 
are unburned (i.e. depleted in Fe-group elements and also IME to 
avoid Ca and Si high-velocity absorption), somewhat richer in IME 
(Si mass fraction 8% in this model) and finally in IME/Fe-group 



elements (Si mass fraction: 15%, Fe-group mass fraction: 9%), re- 
spectively. C II A6580 appears as a prominent line in the spectra. It 
is reproduced using surprisingly low mass fractions of C (3% in the 
outermost zone, 2% everywhere else), as the high luminosity of the 
SN leads to an ionisation state favouring C II over C I. 

Our synthetic spectra are generally too red. A bluer colour 
would have implied a high Tr throughout the atmosphere. We can- 
not implement this, as a low enough Tr is crucial for a reasonable 
ionisation balance, which is needed to match e.g. the ratio between 
the Si III A4563 and Si II A6355 features. 

Within the first six days covered by the observations, the spec- 
trum of SN 2009dc does not evolve very much. For both spectra, 
the photospheric velocities are low (9050 km s" 1 at — 9.4 d and 
8300 km s _1 at —3.7 d) compared to normal SNe la before maxi- 
mum. The most prominent lines due to Si II and S II are somewhat 
stronger at —3.7 d, but are still formed in the same atmospheric 
layers as before. This means that a reduction of IME in the outer 
layers has not only an effect on the —9.4 d spectrum, but also on 
later ones. 



3.2.2 +4.5 d model 

The spectrum and the model have now changed considerably. Lines 
produced by IME still gain more strength, as temperatures drop in 
the outer layers of the atmosphere. The cores of the lines of singly- 
ionised IME (such as Si II) remain at relatively high velocities (as 
at the early epochs). Only the low-velocity wings of the absorp- 
tions are caused by the lower layers of the ejecta. In Si II A6355 
and in the S II W-feature (observed around ~ 5400 A), these low- 
velocity wings can unfortunately not be reproduced, regardless of 
the abundances assumed. This may indicate a too strong ionisation 
in the lower layers of the model, or a too high photospheric veloc- 
ity, which however is required in order not to obtain a hot lower 
boundary spectrum. 

The C II A6580 feature in the observed spectrum has shifted 
to much lower velocities by this epoch. This clearly indicates the 
presence of C in the lower layers. The line's strength at high veloc- 
ities, at the same time, is automatically reduced because of lower 
densities and temperatures (the line originates from a highly ex- 
cited level). Another feature indicating a drop in temperature is O I 
A7773, which becomes deeper in the observations and in the model. 
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Figure 3. Sequence of spectral models (red graphs) for SN 2009dc, based on the AWD3det density profile, and observations (black graphs). Identifications 
given for the most prominent features. 
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3.2.3 +12.5 d model 

The +12.5 d spectrum shows a considerable drop in UV flux. Pho- 
tospheric temperatures are lower, but the line blocking in the UV 
also increases. The Fe-group abundances at the photosphere are 
higher than for the previous spectra. Especially Co II and Cr II are 
now effective in suppressing the flux around ~ 3500 A: the ionisa- 
tion degree is lower than at the previous epoch, favouring recombi- 
nation of Co III and Cr III into the singly-ionised state. Despite the 
decreasing ionisation, the synthetic Fe III feature at ~4300 A is 
deep enough to reproduce the observed line, as significant amounts 
of 56 Ni have decayed into 56 Fe. No additional 54 Fe is assumed to 
be present in the models. 

The observed features are all in all reproduced, with some ex- 
ceptions: at ~4500 A, there is a marked absorption and re-emission 
feature in the observations, which is reproduced only at later epochs 
in our models. The other exception is Si II A5972, which is too 
shallow. We suggest that Na I D starts influencing the A5972 fea- 
ture around this epoch. Na is not included in our models, as past 
attempts to fit the feat ure have suffered fro m deviations in the ioni- 
sation of this element jMazzali et all 19971) . 

The Si II A6355 velocity shows a mismatch between models 
and data, which increases with epoch. The observed line moves to- 
wards the red (lower velocities) as time progresses, while the centre 
of the synthetic line does not change that much. We have already 
removed Si from the outermost layers as far as the earliest spectra 
permit. Part of the velocity drop in the observed spectra may actu- 
ally be caused by a feature forming bluewards, which gets stronger 
with time and emits into the Si II A6355 region. Judging from our 
line list, this is probably a Fe II feature, but we do not reproduce it 
with our models. 

3.2.4 +22.6 d model 

At 22.6 d past B maximum, the photosphere recedes into a zone 
rich in Fe-group elements. The most prominent features in the spec- 
trum are due to singly-ionised Fe-group elements with very strong 
lines (marked in Fig.O. The balance of these lines is mostly a mat- 
ter of temperature and density structure. 

3.2.5 +36.4 d model 

Remarkably, SN 2009dc still shows a photospheric spectrum long 
after maximum, with a v p h of 4700 km s _1 at +36.4 d. The quality 
of the spectral fits is relatively good given the late epoch. The slow 
evolution into the nebular phase indicates higher densities in the 
intermediate and inner part compared to normal SNe la. 

The photospheric velocity at +36.4 d is much lower than at 
+22.6 d. This jump is, however, probably not physically signifi- 
cant, and rather an artefact of the photospheric approximation for 
such late epochs: an unusually low backscattering rate in the latest- 
epoch model shows that this approximation begins to break down 
here. 

3.2.6 Abundances 

The abundance stratification of the AWD3det-based ejecta models 
is shown in Fig. [4] 

The outer zone (v > 9650 km s -1 ) only contains small 
amounts of IME and Fe-group elements. This corresponds to the 
weak line blocking and the relatively shallow features in the early- 
time spectra. Due to the relatively high ionisation, carbon mass 
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Figure 4. Abundance Tomography of SN 2009dc based on the AWD3det 
density profile. The " 56 Nin" and "Fen" lines represent the 56 Ni and sta- 
ble Fe abundances directly after explosion, respectively (cf. Appendix I51 
Table iBlY Panel a): abundances in mass space. Panel b): same, in veloc- 
ity space. In panel c ), th e abundances of t he SN la explosion model W7 
llwamoto et afll 19991 cf. lStehle et aIH 2005) for normal SNe la are plotted 
for comparison. Whe n assuming a zero-metallicity progenitor [model W70, 
llwamoto et alj fl999k the most significant change is the absence of 'Ten" 
in the outer layers, represented by the dotted part of the respective curve. 



fractions of some per cent suffice to fit the observed C II A6580 
line, although it is much stronger than in other SNe la. C must be 
present down to ~ 7500 km s _1 in order to explain the observed 
low velocity in the post-maximum A6580 feature. The majority of 
the material in the weakly burned zones, however, consists of O. 

Between ~9000 and 7000 km s _1 IME are highly abundant. 
The mass fraction of Fe-group elements in this layer, between 12% 
and 58%, is constrained by line blocking of the UV flux. 

Although our model sequences extend to late epochs, we 
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mostly explore regions of incomplete burning. At velocities 
< 7500kms~ 1 , finally, Fe-group elements start to dominate. 
These zones are sampled by the +22.6 d and +36.4 d spec- 
tra. Due to degeneracy, the data can be fitted with different Fe- 
group abundances. We choose rather high mass fractions, consis- 
tent with the fact that we see some Fe-group elements already in 
the —9.4 d spectrum. The mass fractions of " 56 Nio" and "Feo" 
(abundances of 56 Ni and stable Fe directly after explosion, respec- 
tively) together constitute 93% by mass at the +36.4 d photosphere. 
"Feo" is only present deep inside our model ejecta, as our mix 
of Fe-group elements has been set up (cf. Sec. |2.2.5t to resem- 
ble the low-metallicit y W7-based nucleosynthesis calculation W70 
diwamotoetalJl999h . 

With respect to the W7AV70 nucleosynthesis (see Fig.|4j;), all 
burning products are found at much lower velocities in SN 2009dc. 
As a further difference, the "abundance mixing" is bit stronger than 
that of W7, as has also be en found in tomography studies of other 
SNe la dStehle et alj|2005t) . In SN 2009dc, we see C in fairly deep 
layers, and need some Fe-group material in the outer zones to re- 
produce the spectra early on. IME material is present in the un- 
burned and Fe-group layers. Except for the outer layers, the sen- 
sitivity of the spectra on the IME mass fraction is limited. Thus, 
the spectra can also be reproduced with some more down-mixing 
of O, at the expense of Si. Stricter limits on the O abundances in 
the lower layers may be deduced in further studies modelling the 
nebular spectra of SN 2009dc. 

The cumulative chemical yields (total mass of C/O, IME and 
Fe-group elements) of all our models are given and discussed in 
Section |6~2l 

3.3 Abundance profiles and assumptions on 
metallicity/reddening 

We repeat the AWD3det-based tomography for a lower reddening 
E(B — V) ~ 0.06 mag and solar progenitor metallicity, respec- 
tively, in order to examine how our results change with these as- 
sumptions. 

When using a lower reddening, the models have a lower lu- 
minosity and the radiation field in their interior is less intense, de- 
creasing the ionisation. This holds even though we use a shorter 
optimum rise time (Sec. 13. U of only 28 d (which tends to increase 
radiation field densities, cf. Sec. 12.2.3b . The models contain more 
singly-ionised material, which is effective in producing lines (IME) 
and blocking the UV flux (Fe-group). Thus, we reproduce the spec- 
tra with somewhat lower abundances of IME in the outer layers of 
the ejecta, and lower abundances of Fe-group elements. Qualita- 
tively, however, the abundance stratification in velocity space re- 
mains the same (Fig. |5ji). The quality of the models is a bit better 
than with the high standard reddening, which tends to exacerbate 
the models' flux excess in the red. 

A high (solar) metallicity, on the other hand, implies that there 
are lo wer limits to the abund ances of Fe-group elements in the outer 
ejecta dAsplund et ai]|2009h . Furthermore, there is a change in the 
typical nucleosynthesis patterns: signifi cant amounts of 54 F e ap- 
pear in layers with incomplete burning i ll wamoto et aTll 19991) . We 
reflect this, as far as possible, in the abundances of Fe-group ele- 
ments in our models (Fig.|5p). The resulting spectra have a slightly 
worse quality: the outer layers with solar Fe-group abundances 
block flux, heating up the atmosphere and increasing the flux ex- 
cess in the red, where the radiation can finally escape. Photospheric 
velocities are quite a bit higher in order to keep the backscatter- 
ing reasonable and the atmospheric temperatures low enough to 
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Figure 5. Influence of reddening and metallicity on the AWD3det-based 
tomography (plots in velocity space). Panel a): abundances inferred for 
the lower limit reddening E(B — V) ~ 0.06 mag and low metallicity. 
Panel b): same for E(B — V) ~ 0.17 mag (normal), but solar metallicity 
(i.e. higher metallicity than in our standard model). The abundances of O, 
Si, "Feo" and " 56 Ni " in the standard model [E(B - V) ~ 0.17 mag, low 
metallicity - cf. Fig. [4] are plotted as thick, lighter lines for comparison. 



form the spectral lines. This shifts the borders of some abundance 
zones. Still, the abundance stratification we infer is again qualita- 
tively similar to the standard, low-metallicity case. 

All in all, we see that our results are somewhat sensitive to 
reddening and metallicity, but the respective changes are not large 
enough to compromise our conclusions. 



4 TOMOGRAPHY WITH A CORE-COLLAPSE-LIKE, 
EMPIRICALLY MOTIVATED DENSITY PROFILE 

In this part, we conduct a tomography with an "empirically mo- 
tivated" density profile ("09dc-exp"). This profile is constructed 
(Sections |4.1||4.2t from two exponentials (i.e. it is a piecewise ex- 
ponential curve) with a flatter slope in the inner part. The outer part 
is inferred from the evolution of line velocities in the spectra, as- 
suming that lines always form at a density roughly constant with 
time (see Appendix lAl. 

Although we do not assume any explosion model in the be- 
ginning, the profile resembles density profiles found in simulations 
of core-coll apse SNe with C-0 star progenitors (e.g. llwamoto et all 
Il994ll2000f) . With an ejecta mass of ~3M , our models based on 
09dc-exp can be seen as an approach to interpret SN 2009dc as a 
core-collapse explosion. 
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Figure 6. Examples of the density profiles for the outer ejecta constructed 
in Sec. I4.ll / Appendix fA] and the W7 profile as reference. All profiles are 
plotted for a reference time to = 100 s after the explosion, although they 
have been created assuming different rise times t r , 



4.1 Empirically inferred outer density profile and optimum 
rise time 

4.1.1 Grid of empirically-inferred density profiles 

We have constructed a two-parameter grid of exponential density 
profiles for the outer layers of the ejecta (Appendix |Aj examples 
see Fig.[6]l. One parameter is the absolute density scale, which we 
parametrise here in terms of the velocity v x at which the profile 
first intersects the W7 model (cf. Fig. [6]l. For this we allow values 
of 9000, 10000, 12000, or 14000 km s" 1 . The other parameter is 
the rise time (21 d, 22.5 d, 25 d, 27.5 d, 30 d, 35 d or 40 d), with 
which the slope of the density profile inferred from the line velocity 
evolution changes somewhat, but not drastically (cf. Appendix [At, 
The set of exponential profiles, which are given names reflecting 
v x (e.g. "exp9" matches W7 at v x = 9000 km s _1 ), covers well 
the physical parameter space in which solutions for SN 2009dc will 
be found (see below). 

For each v x , we infer an optimum rise time separately. In the 
end, we choose one combination of (v x \t r ), which promises opti- 
mum results, in order to construct the full 09dc-exp density profile. 



4.1.2 Density/rise time test 

The optimum rise time for each v x is inferred from the position 
mismatch of the Si II A6355 feature, comparing the observed spec- 
tra to the synthetic ones. This procedure has already been described 
in Section [3~T| Again, the influence of solar metallicity, and of a 
lower-than-standard reddening is explored by repeating the mod- 
elling for these cases. 

Table [2] gives - for each assumed reddening, metallicity and 
v x separately - the measured line position mismatches between 
the observations and the four optimum models (for different t r ), 
as well as the resulting t r , opl . In Figure [7] we plot spectra for a se- 
lection of these models (for E(B — V) ~ 0.17 mag and negligible 
metallicity). 

Using a smaller v x for the empirical density leads to smaller 
densities throughout the envelope. Looking from the outer layers 
inwards, high densities are reached later, and the zone favoured for 
Si II line formation tends to be deeper inside (which reduces line 
velocities). Thus, the t, values needed to reproduce the low line 
velocities of SN 2009dc are usually shorter for a low v x . If v x 
is very low, the zone dense enough for line formation may be far 
inside where the radiation field is strong. This then favours Si III 



Table 2. Rise time estimates from our spectral fits, using different expo- 
nential density profiles for the outer ejecta (see text), for different values of 
the reddening E(B — V) and metallicity. The table is analogous to Table 
[T] with the difference that several density profiles are tested, corresponding 
to different v x (the slope of the density also changes with i r , cf. Appendix 
lAl however only slightly). The optimum rise time 4 r ,opt is again calculated 
from a linearly-fitted AA — t r relation. 
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35.0 
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with the result that the Si II A6355 line cannot be well formed. For 
this reason, the exp9 models allow us to reproduce the observed 
Si II line d epths only with unre alistically high Si mass fractions 
> 70% (cf. llwamoto etal.lll999l) . if at all (see t r = 25.0 d model 
in Fig.|7](. We therefore expect the actual densities in the outer lay- 
ers of SN 2009dc to be higher than those in the exp9 models. The 
expl4 models, on the other hand, need a long rise time to reproduce 
the Si II A6355 line (TableO and have an excessively high opacity. 
Thus, they constitute an upper limit to the real densities. 
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Figure 7. Models for the — 9.4 d spectrum of SN 2009dc based on the empirically-inferred density profiles [analogous to Fig. [2] for negligible progenitor 
metallicity and E(B — V) = 0.17 mag] . Each of the four panels contains two models for the same density parameter v x and the most extreme t r values listed 
in Table[2] An optimum model would fit the velocity of the Si II A6355 line (marks; inset). With a very short rise time, the velocity (blueshift) of the synthetic 
feature is too high (blue graph); with a very long rise time, the velocity is too low (orange graph). For the exp9-based model with tt = 25 d, the synthetic Si II 
A6355 line lacks depth due to a low density. 



With the intermediate v x = 10000 - 12000 km s~~\ we ob- 
tain the most reasonable models for the —9.4 d spectrum. 



4.2 Construction of the complete 09dc-exp profile 

Therefore, we now conduct a tomography experiment with a "09dc- 
exp" density profile, whose outer part has a v x of 11000 km s _1 . 
As a rise time we use 29 d, which is the mean of the optimum values 
for the explO and expl2 profiles (28 d and 30 d, respectively, see 
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above) under our standard assumptions [E(B — V) = 0.17 mag, 
low progenitor metallicity Z ~ 0]. 

In Fig. [TJ we have plotted the 09dc-exp density. To construct 
this profile, we have used the density profile with v x — 11000 
km s -1 for the outer ejecta. In the inner ejecta, however, we assume 
the slope to be flatter, which is a common feature in SN exp losion 
models (e.g. llwamoto et alJI 1994 1 19991 : IPakmor et alj|201 lh . This 
is done by attaching below some transition velocity an exponential 
inner part with a smaller slope parameter 6j lmcl = 0.1 x 6 0U ter (cf. 
parameter b in Appendix lAl Formula [A2l here 6 out ci = 5.79 x 10 -9 
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significantly with respect to the AWD3det-based models. This will 
be discussed further in Sec. 16. 21 



Figure 9. Abundance stratification of SN 2009dc from our analysis 
based on the 09dc-exp density profile. The abundances of O, Si, "Ferj" 
and " 56 Ni " in the AWD3det-based model [E(B - V) ~ 0.17 mag, low 
metallicity - cf. Fig.|4) are plotted as thick, lighter lines for comparison. 



km -1 s). The transition velocity is chosen such that the total 
mass is 3 Mp, min i mally larger than the ejecta-mass estimate of 
iTaubenberger et al] d201 ll . their Sec. 5.1). The flatter slope in the 
inner zone leads to moderate core de nsities, comparable to those in 
core-collapse explosion models [e.g. IWooslevetal] l ll995h - from 
their Fig. 12 / model 20 we obtain, for 100s past explosion, densi- 
ties of ~ 2 gem -3 at v < 5000 km s -1 ]. 



4.3 Tomography 

4.3.1 Spectral models 

The optimum synthetic spectra calculated with the 09dc-exp profile 
closely resemble those obtained with the AWD3det model (both 
sets of spectra are shown in Fig. [8}. At early epochs (—9.4 d, 
—3.7 d, +4.5 d), the low-velocity wing of Si II A6355 is repro- 
duced better with 09dc-exp. This is due to the lower densities in 
the outer layers (Fig. [TJ, due to which backwarming is less effi- 
cient; temperatures are then lower throughout the atmosphere, so 
that some Si II is present also in the inner parts (where tempera- 
tures are normally highest, so that only Si III exists). 

Later on, the most significant difference between the models 
is the better fit to O I A7773 with 09dc-exp. The AWD3det model 
has too high density at high velocities, causing too much absorption 
in the blue wing of O I A7773. 

In addition, the 09dc-exp model has smaller photospheric ve- 
locities at the latest epochs modelled (+22.6 d and +36.4 d). The 
backscattering rates for these models are therefore no longer in- 
consistent with the assumption of a photosphere, as they were for 
AWD3det. 



4.3.2 Abundances 

The abundance structure inferred with the 09dc-exp profile (Fig. [9} 
is similar to the AWD3det standard case, except for some shifts in 
the zone borders. Photospheric velocities are insignificantly higher 
in the beginning, which corresponds to the assumption of a smaller 
tr. As mentioned above, the opposite holds after B maximum. 

Due to the different density profile, the integrated yields (total 
masses of unburned material, IME and Fe-group elements) change 



TOMOGRAPHY ASSUMING AN EXTERNAL 
CONTRIBUTION TO THE LUMINOSITY 



As already discussed by ITaubenberger et alj J201 ll) . an explosion 
of a single Mch,non-rot WD is most probably inconsistent with 
SN 2009dc if all the observed light is generated internally by Ni 
decay. Such an explosion may, however, explain the SN if strong 
interaction of the ejecta with some circumstellar medium (CSM) 
generates some of the light received from the SN site. In this case, 
the SN needs to contain only a smaller amount of 56 Ni, and the long 
rise time may not be the result of a large progenitor mass. Instead, 
it may be due to an increased optical depth caused by the swept up 
CSM, or due to a longer time-scale on which the interaction pro- 
duces light. 

In an approach inspired by lHamuv et al. I J2003h . we test 
whether the spectra of SN 2009dc can be fitted with less luminous 
models, augmented by some continuous distribution of radiation 
[09dc-int (interaction) model]. For this experiment, we reduce the 
bolometric luminosity in the simulation by 50% at —9.4 d (with re- 
spect to the AWD3det-based model). At the subsequent epochs, we 
keep the absolute luminosity reduction the same (1.5 x 10 43 erg), 
assuming no significant increase of the continuum. From +12.5 d 
on, we begin a slow transition towards normal luminosities, so that 
at +36.4 d the simulation produces 67% of the luminosity in the 
models presented before. The mismatch between the observations 
and the model spectra with reduced luminosity is automatically fit- 
ted with a third-order polynomial after each code run in the fitting 
process. The final model parameters lead to an optimum fit of the 
observed spectrum by the sum of the model and the automatically 
generated polynomial. 

5.1 Density profile and rise time 

The density profile used for these calculations is "w7e0.7" (Fig.|TJ, 
a rescaled W7 model with Chandrasekhar mass, but a kinetic en- 
ergy reduced to 70% of the original value. The smaller kinetic en- 
ergy favours line formation at low velocities and may be interpreted 
in terms of a deceleration of the ejecta by the CSM. It is technically 
implemented by rescaling each point in the W7 velocity-density 
grid according to: 



P 



«W7 



PW7 ' 



K 

-Ek,W7 

K 

■Ek,w7 



1/2 



-3/2 



0.7. 



where (t>w7, Pwi) are the W7 values and r Bk 

The interaction-based model is more of an exploratory nature: 
as the scenario implies a deceleration of the ejecta when interacting 
with the CSM, the assumption of force-free expansion is not really 
accurate here; the force-free time evolution of the density assumed 
by the spectrum code will deviate to some degree from that of real 



4 The difference between the observed and the synthetic spectrum is fitted 
using a Levenberg-Marquardt method. The wavelength range for the fit is 
restricted to < 7500 A so that the procedure is consistent among the spectra 
whose IR extent varies. 
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Figure 8. Sequence of spectral models for SN 2009dc, based on the 09dc-exp density profile (red graphs). The models based on the AWD3det density profile 
(cf. Fig. [5) are plotted in blue for comparison, as well as the observations of SN 2009dc (black graphs). Identifications are given for prominent features. 
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interacting ejecta. We therefore want to draw only qualitative con- 
clusions from this part of the study. 

In line with this, we have refrained from inferring the rise time 
in such detail as in the other cases. Instead, we have taken advan- 
tage of the fact that the outer part of the w7e0.7 profile is quite sim- 
ilar to the 09dc-exp profile (cf. Fig.[T](, for which we had inferred a 
rise time of 29 d (Sec. 14. lb . Starting with this value, we began mod- 
elling the earliest spectrum and obtained a somewhat too red Si II 
A6355 line. We then reduced the risetime of our model in crude 
steps of 1 d. We finally obtain a good match of the observed and 
synthetic line for a rise time of 26 d. 

5.2 Spectral models 

Our synthetic spectra for the 09dc-int model, with their respective 
intrinsic and continuum contributions, are shown in Fig. [TO] The 
fit to the observations is very good, especially at the early epochs 
where the continuum contribution now allows us to perfectly match 
the blue colour of the observed spectra. The polynomial continua 
become less luminous and redder with time. 

The most marked mismatch between synthetic and observed 
spectra remains in the ~ 5700 A feature at later epochs, where we 
have not attempted to fit the Na I D line. Furthermore, the Si II 
A6355 line keeps too strong a blue wing at late epochs, even more 
than in the AWD3det and 09dc-exp models. This obviously de- 
pends on densities and ionisation in the models, and it is - as other, 
minor mismatches - not a major concern at the current stage. 

We emphasise that our continua (as we have tested) can not 
well be fitted by black-body curves. We speculate that the shape 
may be explained by emission of overlapping lines e.g. of Fe- 
group elements. In order to test whether this is realistic, we tried 
to replace the polynomial continuum by a spectrum of the Type 
Ibn SN 2006jc [taken about three weeks past maximum light, 
IPastorello etall <2007h l. We cut out the He emission features from 
this spectrunjf]. What remains is a pseudo-continuous flux distri- 
bution, which at the epoch chosen should be - at least in part - 
due to atomic lines excited by interac t ion of the SN ejec ta with a 
CSM frolev et al.|[2007l ; IChugaill2009l ; I Smith et al. 1 120091) . Indeed, 
the spectra of SN 2009dc at -9.4 d, -3.7 d and +4.5 d are well 
reproduced by adding the "modified" SN 2006jc spectrum (multi- 
plied by a scaling factor which is the same for all three epochs) to 
our photospheric model (Fig.[l0] upper three panels - dark blue and 
brown graphs). At +22.6 d and +36.4 d, our continua are some- 
what redder than any of the spectra of SN 2006jc presented by 
IPastorello et alj d2007l) . The pseudo-continuum in SN 2009dc ob- 
viously is generated by less and less energetic processes with time. 
To understand what information this gives on the ionisation and 
excitation state of the medium, and on the interaction process go- 
ing on, more extensive radiative-hydrodynamics modelling will be 
needed. 

5.3 Abundances 

Despite the physical differences, the abundance structure (Fig. II It 
remains similar to the one obtained with the AWD3det model 
(again in velocity space - not in terms of the integrated yields dis- 
cussed later in Sec. 16. 2b . The increase of burning products from the 

5 By this we mean inserting a constant flux at 4375-4550 A, 
4900-5125 A, 5750-6050A, 6650 - 6700A, 7000-7150A and 
7200-7350A. 



100 



10 




0.01 1 ' ' • — : 1 

5000 10000 15000 

velocity [km/s] 



Figure 11. Abundance Tomography of SN 2009dc assuming 
Chandrasekhar-mass ejecta, with some of the light generated by in- 
teraction with a CSM: ()9dc-int model (abundances plotted in velocity 
space). The abundances of O, Si, "Feo" and " 56 Nio" in the AWD3det- 
based model [E(B — V) ~ 0.17 mag, low metallicity - cf. Fig. |4) are 
plotted as thick, lighter lines for comparison. A metallicity of 0.5 X Z@ 
was assumed when creating the model (cf. Sec. 13.31 . 

outside inwards is a bit steeper. Higher abundances are needed to 
fit the observed line depths as the continuum contribution makes 
the lines shallower in the final model (sum spectrum). We have 
assumed an intermediate metallicity here (half-solar) because in 
09dc-int, according to test models we calculated, the reproduction 
of the spectra does not depend much on the metallicity assumed. 



6 DISCUSSION 

With our tomography models, we tested different possible explo- 
sion scenarios for SN 2009dc. Here, we comment on the viability of 
these, summarising the constraints under which they are consistent 
with our results. We separately discuss fully-radioactive (Sec. (6/3} 
and interaction-luminosity (Sec. |6.4t models. As a starting point, 
we discuss the limits on the rise time (Sec. l6.lt and on the chemi- 
cal yields of SN 2009dc we obtain from our spectral models (Sec. 
16721 . 

6.1 Lower limit on the rise time of SN 2009dc 

In Sections [3.1ll4.1l and l5.ll we have inferred best-fit rise times for 
SN 2009dc assuming different density models, which well cover 
the physically sensible parameter space for the SN ejecta. Without 
assuming interaction, the optimum rise time inferred is lowest with 
the exp9 or explO profiles; the respective values constitute lower 
limits to the real rise time of the SN for this case. For a reddening of 
E(B — V) ~ 0.17 mag and negligible progenitor metallicity, we 
derive (from Table |2): 

tr, 09dc > 25 d 

Here, a possible error of 3 d (cf. Sec. 13. lb has again been ac- 
counted for already. The lower limit also holds for solar progeni- 
tor metallicity, as the respective models need somewhat longer rise 
times. With a weaker extinction of E(B — V) ~ 0.06 mag, the rise 
time limit inferred from the models is t r ,09dc > 21 d (but, as already 
mentioned in Sec. 13. II there is no compelling reason for such an 
assumption). 
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Figure 10. 09dc-int models for SN 2009dc. Each observed spectrum (black) is modelled as the sum (orange) of an intrinsic SN model spectrum (red) and a 
continuum (light blue), resulting from interaction with circumstellar material. The continua (see also text) are third order polynomials which are used only 
where their values are greater than zero and bluewards from their minimum (which occurs at > 7000A). For the earliest thre e epochs, the continua are very 
similar in shape to a post-maximum spectrum of the interaction-dominated SN 2006jc with the He emission lines removed i Pastorello et al. 2007t see also 
text). Adding the SN 2006jc spectrum (scaled; dark blue lines) as a pseudo-continuum to our photospheric models therefore reproduces the observations very 
well (brown lines). 
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Table 3. Integrated yields of C+0, IME and Fe-group elements in the mod- 
els. The Fe-group values include the respective amount of 56 Ni. For the core 
below the +36.4 d photosphere [mass M (core)], we cannot constrain the 
abundances. We expect it to contain a significant amount of 56 Ni. 



model 


M(C+0) M(IME) A/ (Fe-group) M( B6 Ni) M(core) 




[M ] 


[M ] 


[M ] 


[M ] 


[Mgl 


09dc-AWD3det 


0.92 


0.28 


0.30 


0.39 


0.30 


09dc-AWD3det, low E(B — V) 


1.1(1 


0.34 


0.3-1 


(1.31 


0.23 


09dc-AWD3det, solar metallicity 


0.75 


(1.33 


0.62 


(1.46 


0.30 


09dc-exp 


0.57 


(1.63 


1.38 


1.21 


(1.39 


09dc-irjt 


0.38 


0.3 1 


0.36 


(1.31 


(1.33 



In the same manner, one can read an upper rise time limit of 
~ 36 d from Tabled; this limit is however less stringent than the 
respective observational limit [30 d. ITaubenberger etal I fcOllh l. 

In the interaction scenario, we found the the optimum rise time 
to be 26 d (Sec. 15. lb . i.e. within the given range. 

6.2 Limits on the integrated nucleosynthesis yields 

The total masses of C+O, IME, Fe-group elements and 56 Ni we 
find in our different tomography models are compiled in Table [5] 
The 56 Ni masses shown do not include the 56 Ni in the opaque 
core of the ejecta below the +36.4 d photosphere. The mass of this 
core, which is expected to contain a significant additional amount 
of 56 Ni, is also given in the Table. 

6.2.1 Mass of 56 Ni: Arnett's rule vs. spectral analysis 

The 56 Ni mass allowed by each tomography model is an important 
indicator for the consis tency with SN 2009dc, as we can compare 
it to what Arnett's rule dArnettll 19821) requires to produce the light- 
curve maximum of SN 2009dc. This rule states that at maximum 
light the luminosity roughly equals the instantaneous energy input 
from radioactive decay: 

£(tmax) = -MNi-56(imax) X C^Ni-56 + Mco-56 (£max) X £q>56 (1) 

(x is the energy gained from radioactive decay of an isotope 
X per unit mass (of this isotope) and unit time. The time- 
dependent masses of 56 Ni and 56 Co are calculated using the half- 
lives ti/2.Ni-56 = 6 . 077 d and ti/2,co-56 = 77.27 d. Applied as in 
ITaubenberger et alj fcOllh . the rule yields a 56 Ni mass of 2.0 ± 0.7 
or 2.3 ± 0.8 Mq for an assumed rise time of 25.0 d (our lower 
limit) or 30.0 d (the observational upper limit), respectively. The 
error estimates include uncertain ties in the distance, extinc tion and 
in applying Arnett's model (cf. ITaubenberger et al]|201 ll) . These 
values for the 56 Ni mass, which are significantly larger than the 
Chandrasekhar mass Mch.non-roi, may be inaccurate in case of large 
asym metries in the supernova dHillebrandt et al.ll2007l ; Isim et all 
120071) . As spectropolarimetry of SN 2009dc however do es not in- 
dicate extreme global asymmetries dTanaka et al.ll2O10l) . one can 
however say that Arnett's rule in combination with our rise time 
limit practically rules out a Chandrasekhar-mass progenitor if 56 Ni 
is the only source of light in SN 2009dc. 

The amount of 56 Ni in our tomography models, which we can 
compare with the values given, is somewhat uncertain due to degen- 
eracies especially at late epochs. We have chosen 56 Ni abundances 
as high as possible. If we additionally assume the opaque core of 
the ejecta at +36.4 d to be dominated by 56 Ni, we can perform a 
consistency assessment which is conservative in the sense that we 
will not reject models when uncertain. 



The smallest 56 Ni mass (< 0.7 Mg) is derived for the 09dc- 
int model. This model is a special case as it assumes part of the 
luminosity not to come from 56 Ni. Comparing to the 56 Ni mass re- 
quired by Arnett's rule, one recognises that the interaction with the 
CSM would have to produce ~ 2/3 of the maximum luminosity. 
This is somewhat more than we assumed in the models, but taking 
into account uncertainties the model will be more or less consistent. 

For models assuming the luminosity to be intrinsic, we obtain 
sharper constraints. The AWD3det based models only allow for a 
56 Ni mass smaller than 0.8 Mq, which is a severe problem if one 
wants to explain SN 2009dc. We further discuss this in Sec. 16. 3. 1~| 
The 09dc-exp model has a larger mass concentrated in the central, 
56 Ni-rich layers. Therefore it contains up to 1.6 Mq, which is con- 
sistent with Arnett's rule within the error bars. 

6.2.2 Minimum mass of directly-synthesised Fe, IME and C+O 
inSN2009dc 

Our models explore a range of explosion scenarios expected to ex- 
plain SN 2009dc, which imply different density structures and dif- 
ferent conditions for ionisation and excitation in the ejecta. We can 
therefore regard the minimum of the directly-synthesised Fe ( 56 Fe 
from 56 Ni decay excluded), IME and C+O masses found in our 
models as rough lower limits to the actual nucleosynthesis yields 
oftheSN. 

The amount of directly-synthesised Fe (which is roughly the 
difference between the Fe-group and the 56 Ni values in Tabled 
is practically insignificant (mostly smaller than 0. 1M©). For the 
intermediate-mass elements, we obtain 

MlME,09dc > 0.3 M , 

which is the IME mass in the 09dc-AWD3det model. 

For the C+O mass, to which C contributes only little, we ob- 
tain from the 09dc-int model 

M c +o,09dc > 0.4 M . 

The oxygen included here is needed for reproducing the observed 
O I A7773 feature. 

6.3 Candidate explosion models without CSM interaction 

6.3.1 Rapidly-rotating WD models 

With our AWD3det-based models, we explored the possibility of 
explaining SN 2009dc as an explosion of a massive, rotating C-O 
wh ite dwarf. Rotating w hite dwarfs with up to ~ 2 M© may ex- 
ist dYoon & Langerll2005h . although there is no observational proof 
and the m ass limit is a matter o f debate [see the sum marizing dis- 
cussion of IScalzo et ail feOlOh l. lHachisu et al.1 J20I2I) propose that 
even masses up to ~2.5 M© may be realistic. 

Assuming the AWD3det density, our spectral models do not 
allow enough 56 Ni in the intermediate and outer (i.e. high-velocity) 
zones of the ejecta in order to obtain the large 56 Ni masses expected 
for SN 2009dc (cf. Sec. 16.2. It . Nucleosynthesis calculations for 
the hydrodynamical explosion model indicate a 56 Ni mass of ~ 
1.5 M©, which would be compatible with SN 2009dc. However, 
such efficient nucleosynthesis within the white dwarf implies (in 
terms of energetics) that burning products, in particular Fe-group 
elements, are ejected with high velocities, contradicting our models 
(Fig.|T2]l. Also, the abundances of IME and C/O in the AWD3det 
simulation (< 0. 1 Mq, respectively) are lower than the limits given 
in Sec. 16.21 
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Figure 12. Comparison of the abundance stratification inferred from the 
spectra (using t he AWD3det d ensity) with t he nu cleosynthesis of the 
AWD3det model iPfannes et all feOlObl) ; iFinkl feOld) . We show the total 
abundances of three element groups (C+O, IME and FeG / Fe-group ele- 
ments), as they are used in the AWD3det calculations. The AWD3det mass 
fractions have been obtained from the original 3-D model by adding up the 
masses of C+O, IME and Fe-group elements, respectively, into radial ve- 
locity bins. The lower limit of the velocity axis corresponds to the +36.4 d 
photosphere. 



some time, and the primary will slowly be burned to an O-Ne-Mg 
WD. This WD will subsequently undergo accretio n-induced col- 
lapse jNomoto & Ibenll 19851 : Iskio & Nomotoll 19981) . Therefore, in 
order to explain a SN explosion with a non-violent merger, C igni- 
tion has to be avoided. This may be possible at least in a few con- 
figurations l Yoon et al J2007t) . The primary may then grow towards 
its Chandrasekhar mass and explode. 

In any case, if a "merger SN la" is made by two WDs with 
significantly different mass, the 56 Ni mass generated will not be 
larger than Mch,non-rot, as would be needed to explain the light curve 
of SN 2009dc: the material of the secondary will b e too dilute to 
be co nverted into 56 Ni when the explosion occurs (cf. lPakmor et al.l 
l2012h . WD pairs of roughly equal mass and a total mass ~2Mq, in 
turn, are probably not too frequent (the exact occurrence rates are 
not exactly known). 

Therefore, we suggest that WD mergers will rather not pro- 
duce an extreme 56 Ni mass, as it is required for SN 2009dc 
without assuming CSM interaction. It seems more likely that a 
double-degenerate system could explain SN 2009dc if an extended, 
C/O-rich circumstellar medium was generated during the accre- 
tion/merging process, and some of the luminosity of SN 2009dc 
was contributed by an interaction of the ejecta with that medium 
(see Sec.|6.4|l. 



The major uncertainty in our analysis - besides the maximum 
WD mass - is that we assume spherical symmetry in our spec- 
tral analysis and in applying Arnett's r ule. The AWD3det eject a 
are significantly asymmetric dFinkll20ld: cf. lPfannes et' al. 2010b), 
possibly provoking viewing-angle effects in the spec trum. Another 
cavea t is that explosions of extremely massive WDs dHachisu et al .1 
120121) may produce the 56 Ni mass for a bright light curve without 
having too-high-velocity Fe-group material. Therefore, we can not 
yet completely exclude SN 2009dc to be explained by a rotating 
WD model, although the factor > 2 between the 56 Ni mass in our 
AWD3det-based model and the 56 Ni mass calculated with Arnett's 
rule (Sec. l6.2Tt seems difficult to overcome. In the end, 3-D radia- 
tive transfer simulations of massive WD explosions may allow for 
a final judgement. 



6.3.2 Mergers of C-0 white dwatfs 

Mergers of two C-0 WDs are another channel for creating a "super- 
Chandrasekhar- mass object" which may undergo a thermonuclear 
explo s ion (e.g. | Yoon et all 2007; Pa kmor et alj|2010al : iDan et al .1 
|2012|) . Ifficken et alj d2007l) and lScalzo et al.l j2010h have already 
considered merger scenarios to explain the SC SNe la 2006gz and 
2007if. 

For small mass ratios (i.e. if one WD is much heavier than 
the other), accretion is a self-stabilising and relatively slow process 
dD'Souza et al ]|2006l : iMotl et al.ll2007l) : the separation of the white 
dwarfs gets larger when mass is accreted, stopping the accretion. 
Such systems may make "normal" Chandrasekhar-mass SNe la if 
the primary (the heavier white dwarf) ever grows sufficiently. 

For higher mass ratios, for which mass transfer is unstable, 
the most probable outcomes are a violent merger or an accretion- 
induced collapse. In a violent merger, the WDs are quickly dis- 
rupted or at least heavily deformed. Du ring this process, th e igni- 
tion of an off-center detonation is likely dPakmor et al.l201ll) . and a 
SN la explosion results. If conditions are not extreme for a violent 
merger (i.e. no detonation is ignited), C will ignite off-centre after 



6.3.3 Core collapse and Type I— SNe 

The 09dc-exp model, which has an ejecta mass of 3 Mq, corre- 
sponds to an explosion of an extended progenitor star instead of 
a WD. This would imp ly either core-collapse, or a thermonuclear 
explosion within a star dErgma & Paczvns"killl974l) . In any case, the 
explosion model has to produce an ejecta mass > 2 Mq and a 56 Ni 
mass > 1.5 Mq (cf. Sec I6.2t 

lUmeda & Nomotol d2008l) discussed the possibility of produc- 
ing 56 Ni masses of some Mq in a core collapse. Their models, 
based on stars with main-sequence masses of M ^ 25 Mq, yield 
Ni masses as we need them for SN 2009dc. However, they still 
contain H and He envelopes (Mh+hc > 10 Mq), which would have 
to be removed by stripping or mass loss in order for the SN to 
potentially display a T ype la spectrum. One cha racteristic feature 
of the (1-D) models of lUmeda & Nomotol d2008h is an abundance 
stratification (Fig. 5 of their paper) which shows similarities to our 
tomography results and generally to SN la models. O is not signifi- 
cantly mixed with 56 Ni on macroscopic scales. This may disfavour 
the appearance of nebular O lines (which are normally present in 
core-collapse SNe), and therefore ensure co mpatibility with the 
non-d etection of these lines in SC SNe la dTaubenberger et all 
l201ll) . It is necessary to explore wh ether envelope stripping w ould 
significantly alter the conclusions of lUmeda & Nom oto (2008)), and 
whether the general structure of the model (abundances, densities) 
would remain the same in a multi-D core-collapse simulation, in 
which th e explosion energy is d eposited in the med ium by neutri- 
nos (e.g. iMarek & Jankall2009l) or acoustic waves dBurrows et al .1 
120071) . 

Thermonuclear explosions within stars, on the other hand, 
have long been proposed but never been verified as a model for 
SNe dErgma & Paczvnskill974l ; llben & Renzinill983l) . Stars with a 
zero-age-main-sequence mass < 8 Mq form degenerate C-O cores 
after core He burning. After mass loss in the asymptotic-giant- 
branch (AGB) phase, the C-0 core usually remains as a bare white 
dwarf. In the most massive AGB stars, however, the degenerate core 
may grow towards the Chandrasekhar limit, and an explosion may 
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commence. The explosion would eject up to a Chandrasekhar mass 
of 56 Ni, corresponding to the entire stellar core. The composition 
of the ejecta depends on the progenitor material and its nuclear re- 
processing. In case only Fe-group material is produced in the core, 
a significant amount of C, O and IME wou ld have to be produced 
in the envelope (cf. the limits in Sec. l6.2.2] l. llben & Renzinildl983h 
however suggested that H and He may be visible in the spectra 
of such SNe, and therefore dubbed them "Type IV2" explosions 
(i.e. thermonuclear with a core-collapse-like display). This must 
be avoided, and it is not clear whether this is realistically possible 
e.g. by nuclear reprocessing of the envelope. Small quantities espe- 
cially of H e may also be present without lines being detectable in 
the spectra dHachinger et al.l2012l) . In the end, the Type F/2 scenario 
has to be simulated with recent hydrodynamic or radiation-hydro 
codes in order to make more definite predictions on its nucleosyn- 
thesis and observational outcome. 



6.4 Candidate explosion models with CSM interaction 

Finally, we discuss a scenario which involves light generation by 
significant interaction of the ejecta with a CSM. The 09dc-int 
model (Sec.|5j shows that the spectra can successfully be explained 
with light from 56 Ni plus an "external" light contribution. 

The external contribution in our models is continuous, but 
does not have a black-body shape. It may be explained by emit- 
ting lines which overlap: at the early epochs the pseudo-continuum 
is very similar t o a sp ectrum of the interacting SN Ibn 2006jc 
dPastorello et al.l l2007t) some weeks past maximum light, with 
its He emission cut out. Such SN Ibn spectra have been sug- 
gested to have a pseudo- continuum from overlapping e mission 
lines of Fe-group elements dFolev et al .l2007l.lChugail2009l : cf. also 
ISmith et al.ll2009h . 

Presumably, a significant optical depth of the emitting zones is 
needed to allow for "partial thermalisation" and generate a smooth 
flux distribution without major glitches. For the light of the SN to 
pass to the observer in such a scenario (as indicated by the relatively 
normal SN spectra), the emitting material has to be concentrated, 
e.g. in clumps or in a torus, so that a major fraction of the lines of 
sight between SN and observer remain free. 

With the additional light contribution, one may explain "SC" 
objects as an explosion of a WD with a mass of roughly Mch,non-rot. 
The CSM itself has probably to consist of C, O, or heavier ele- 
ments, as no H and He lines are visible in the photospheric or 
nebular spectra of SC SNe la. Small en ough quantities of H e 
may again be present without making lines dHaching eretalj2012l) . 
The CSM has to be sufficiently close to be reached by the ejecta 
within a timescale of 10 d (i.e. , for typical ejecta velocities of 
~ 10000 km s" 1 , at a distance of < 10000 km s" 1 x 10 d ~ 10 15 
cm). The light production has to continue for several weeks; if this 
required sweeping up new material continuously, the extent of the 
CSM would have to be at least of the same order of magnitude. 

It is worthwhile to think of progenitor systems possibly gen- 
erating such a CSM. One model to consider are non- violent WD 
mergers with unstable accretion. If the currently predicted scenario 
of burning to an O-Ne-Mg WD and accretion-induced collapse was 
avoided for such a system (cf. Sec. 16.3. 2l . the primary WD could 
explode when reaching Afch,non-mt (or also a higher mass in case 



6 There is only a very short period of mass transfer before a violent merger 
and there is no indication that it could generate an extended CSM (e.g. 



iPakmor et alj2010d,l20nll2012b . 



of fast rotation, for which however tidal locking would have to be 
inefficient). The extent of the CSM, which would have to be pro- 
duced before the explosion, would probably hav e to be much larger 
tha n found in simula tions of such systems by lYoon et al] d2007l) 
and lDanetail d2012h . An alternative possibility would be accre- 
tion from a star which has lost most of its envelope, such that a 
H-deficient and maybe C/O-rich CSM could form. In any case, the 
explosion in an interaction-based scenario has to produce a signif- 
icant, but not an extreme amount of 56 Ni (just as a normal SN la, 
cf. Sec.[5j. 

The CSM interaction scenario has several appealing aspects. 
It may not only explain the luminosity of SC SNe la, but also a 
slow-down of the ejecta and therefore the low line velocities. The 
w7e0.7 density model allowed for good fits to the spectra (09dc-int 
models). In this model, 30% of the kinetic energy of the explosion 
were assumed to be transferred to the circumstellar shells. Crudely 
taking the distribution of kinetic energy in the ejecta to be uniform, 
this would mean that the circumstellar material would constitute 
30% of the total mass (i.e. 0.6 Mq when assuming a Afch,non-mt 
primary). 

Only a small part of the energy lost to the CSM (which in total 
is roughly 0.7 x 10 erg, i.e. 30% of E^wt), would have to be con- 
verted into luminosity with time: the radiated energy differs only by 
some 10 49 erg between normal SNe la and SC objects. Eventually 
the complete CSM would be dragged with the expanding ejecta, 
and the interaction would become less violent. The contribution to 
the luminosity would then be expected to gradually vanish, which 
is consistent with the continuum getting less intense in our models. 

A relatively strong argument against an interaction model is 
the lack of direct evidence in SN 2009dc; in Type Iln and Ibn 
SNe, narrow H/He line emission is observed as a sign of interac- 
tion. Moreover, the light-curve of SN 2009dc after the photospheric 
phase sho ws a slope compatible w ith 56 Co decay for a relatively 
long time dSilvermane7al]|201ll) . suggesting that all light may 
be generated by radioactivity. Whether all this really constitutes a 
problem for an interaction scenario must finally be judged on the 
basis of detailed radiation-hydrodynamics calculations. Also, pop- 
ulation synthesis studies and explosion simulations should clarify 
which progenitor systems for this scenario may exist in significant 
numbers. 



7 SUMMARY AND CONCLUSIONS 

We have conducted an analysis of SN 2009dc, a representative ex- 
ample of the poorly-understood "Super-Chandrasekhar" (SC) SN 
la subclass. These SNe are characterised by an extremely high lu- 
minosity, and by spectral lines typical of SNe la, with normal or 
even low velocities. The photospheric phase lasts unusually long, 
indicating dense ejecta. Using a radiative-transfer code to model 
spectra observed up to one month after B maximum, we have 
inferred the abundance stru cture of the SN ["Abundance Tomog- 
raphy", IStehle et"aL assuming different density profiles. 
Thus we have tested different explosion scenarios for consistency. 

With a density profile repres enting an ex plosion of a massive, 
rotating WD ["AWD3det" model. iFinkfeoifj) and Fink et al. (2012 
in prep.) - an improve d version of the "AWD3 detonation" model of 
IPfannes et ail d2010bh l. the abundance profile we infer allows Fe- 
group elements only at velocities < 9000 km s _1 . As only 1.1 Mq 
of material are ejected at such low velocities, the 56 Ni mass allowed 
with this model is low. In spherical symmetry, as assumed for the 
spectral models, Arnett's rule for the light curve dArnettlll982l) de- 
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mands the 56 Ni mass of the SN to be ~ 2.0 M Q (for t r = 25 d). 
This discrepancy leads us to disfavour a rotating single WD model 
for SN 2009dc. However, the model may be rescued if either the 
explosion is s trongly asymmetric [rendering our initial assump- 
tions wrong -iHillebrandt et 10 J2007h : ISim et al.l d2007h . but see 
iTanaka et al] J2010h l or if explosions of C-0 WDs more massive 
than 2.0 Mq existed dHachisu et alj2012h . for which the ejecta will 
be slower if the amount of burning products is not increased pro- 
portionally to the mass. 

An alternative scenario may be core-collapse or a thermonu- 
clear explosion within a star (Type IV2 SN). This allows for a larger 
ejecta mass, leading to lower ejecta velocities. Thus, a larger 56 Ni 
mass becomes compatible with the spectra. We showed this by con- 
ducting a tomography with an empirically constructed 3 Mq den- 
sity profile. Although the integrated masses of 5f> Ni, IME and C/O 
are different from the AWD3det case, the abundance structure in 
velocity space remains similar. O dominates at v > 9000 km s _ , 
intermediate-mass elements below, and Fe-group elements includ- 
ing 56 Ni, in the core. This stratification would have to be repro- 
duced by a valid explosion model. 

Finally, we explored whether a tomography of SN 2009dc can 
be perfor med assuming a no rmal SN la explosion enshrouded by a 
CSM (cf. lScalzo et al.l 2010h - The interaction of the ejecta with the 
CSM would in this scenario produce additional luminosity, distin- 
guishing "SC" SNe la from normal ones. Thus, the exploding WD 
would not have to produce excessive amounts of 56 Ni, and would 
not need to be supermassive. In order to test this model, we cal- 
culated synthetic spectra with an intrinsic SN luminosity reduced 
by 50 — 33% (with respect to the observed luminosity), using a 
W7 (Afch.non rot) SN la density profile scaled to 70% of its origi- 
nal kinetic energy. At each epoch, the difference between observed 
spectrum and model was fitted with a third-order polynomial, rep- 
resenting the luminosity contribution by interaction. This way, we 
obtained an outstanding reproduction of the observed spectra. The 
polynomial-fit "continuum" we used turned out to be very similar 
to an interaction-dominated spectrum of the Type Ibn SN 2006jc, 
apart from the He emission in that SN. In the context of an inter- 
action scenario, the absence of H and He lines in SN 2009dc can 
only be explained if the companion is H-deficient and He-deficient 
(or not showing He lines because of the ionisation/excitation con- 
ditions). The scenario may thus correspond to a merger of two C-0 
WDs, if both a violent merger and an accretion-induced collapse, 
which currently are considered the most likely outcomes, can be 
avoided. Another possible progenitor system may be a C-O WD 
with a donor star which has lost parts of its envelope. 

We conclude that the popular rotating / supermassive WD sce- 
nario is only one of several candidate models for explaining "SC" 
SNe la. While there are problems with this model, it might work 
with extremely massive WDs or when the explosions are strongly 
asymmetric. The alternatives have been largely unexplored in terms 
of hydrodynamic modelling. Therefore we hope that the present 
work encourages modellers to conduct ab-initio hydrodynamic and 
radiative-transfer simulations of all the scenarios we have dis- 
cussed. The abundance profiles we have inferred may serve as a 
guideline for these simulations. 
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APPENDIX A: EMPIRICALLY MOTIVATED DENSITY 
PROFILES 

We construct part of the density profiles we use in our study in- 
dependent of explosion models (cf. Sec. 14.2b . The idea is to infer 
the outer density profile from the observed time evolution of the 
velocities at which prominent spectral lines form. Such velocities 
are measurable as absorption blueshiftsQ in the spectra. In this ap- 
pendix, we describe in some detail the method to infer the outer 
density, and present the results. 

As an initial step, we make the (very approximate) assumption 
for spectral lines to form at a density constant over time and exploit 
the fact that the ejecta are in force-free, homologous expansion (cf. 
Sec. @r = «xf). With these assumptions, the evolution of line 
velocities is directly related to the radial density profile. The density 
at a space-time point (v2, t2) with respect to the density at another 
point («i,ti) can be calculated from the radial density profile. In 
turn, knowledge on the ratio of the density at two points (v\,ti) 
and (v2,t2) allows us to constrain the shape of the density profile 
[i.e. p(v2,to) can be calculated from p(vi, to) for arbitrary to}. If 
we have velocities of line formation vi...„ at times ti...„, and assume 
a density ratio of 1.0 (constant density) at these points, we obtain: 

p(vi,ti) = p(vi,ti) (constant density) and 
p(vi,ti) — p(vi,ti) x (ti/t\) A (force-free, horn, expansion) 
=> p{vi,t x ) = x {U/hf (Al) 

Here, the time offsets U from explosion onset depend on the as- 
sumed rise time t s . For each spectrum, i, is calculated by adding the 
assumed rise time t Y to the observed epoch relative to B maximum. 
We thus obtain from our measurements a two-parameter family of 
models with t t and the absolute density scale [here appearing in the 
form of p(vi,ti)] as free parameters. 

We perform line velocity measurements in all spectra ob- 
served up to five days past B max. With this choice, we avoid 
the epochs at which the velo cities of some lines become constant 
(cf. iTaubenberger et al]|201 ll) as the photosphere moves below the 
zone where the respective elements are abundant. We use sev- 
eral prominent absorption features due to different species (Si II 
AA4130,6355, Si in A4563, C II A6580). Taking into account vari- 
ous ions with different ionisation and excitation energies and eval- 
uating the average velocity, we make sure that the zone of line for- 
mation we infer does not track the excitation/ionisation conditions, 
but the density (as far as possible). 

The measured velocities are given in Table lAll Exactly speak- 
ing, blueshifts measured in spectra give us velocity components 

7 For the actual measurements, we use the procedure in lHachinger et al] 
i.e. we determine a pseudo-continuum and take the wavelength of 
maximum fractional depth as the blueshifted position. 
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Table Al. Velocities measured in the Si II AA4130,6355, Si III A4563 and 
C II A6580 lines in the observed spectra of SN 2009dc up to 5 d past B 
max. We give the individual velocities, their average and a radial velocity 
of line formation v\j. The values iilf are the V{ used in equation )AU . 



epoch (rel. 
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I'LF 


to B max.) 


(Si II A4130) (S 


i II A6355) (Si III A4563) (C 11 A6580) 


(line avg.) 


(radial) 




[kms" 1 ] [kms -1 ] 


[kms" 1 ] 
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[kms- 1 ] 


[kms- 1 ] 


-9.4 


8298.9 


9071.6 


7061.8 
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8643.6 


10804.5 


-8.5 


7906.2 


9241.8 


6665.3 


10168.7 


8495.5 


10619.4 


-7.5 


7628.1 


8754.1 


6457.1 


9963.1 


8200.6 


10250.8 


-3.7 


7219.0 


8260.8 


5853.6 


9517.6 


7712.7 


9640.9 


2.4/2.5 


6224.3 


7757.0 


5333.7 


7374.8 


6672.5 


8340.6 


4.5 


6021.9 


7482.9 


5197.8 


7191.3 


6473.4 


8091.8 



along the line of sight. The corresponding radial velocity of line 
formation is somewhat higher. Thus, we calculate radial velocities 
of line formation by multiplying the average line velocities with an 
estimated factojfjof 1.25 (see Table lAlt . These radial velocities are 
used as Vi in Equation J A 1 1 > . 

In order to write the inferred density profile (for the outer 
ejecta) in a compact form and to avoid kinks in the profiles, which 
may result from measurement uncertainties, we fit exponentials 

p(v, t) = ax exp(-bv) x (to A) 3 (A2) 

to the measured values p(vi,ti). Here, to is again an arbitrarily 
chosen reference time. As the observations do not constrain the 
absolute density scale (i.e. a in this formula), the fit only gives a 
meaningful b which depends on the assumed rise time t r (t = t r + 
epoch rel. to B max.). For obtaining an actual density profile, t r 
and finally a have to be fixed. We therefore construct a grid of den- 
sity profiles with different t t and a (absolute density scales). We 
assume rise times t, of 21 d, 22.5 d, 25 d, 27.5 d, 30 d, 35 d, and 
40 d and choose a s uch that our models intersect the W7 "stan- 
dard" SN la mode0 jNomoto et all 19841) at 9000, 10000, 12000, 
or 14000 km s _ . As the W7 density is relatively flat in this ve- 
locity range, a smaller intersection velocity v x results in a lower 
a and lower overall densities. This is illustrated in Fig. [6] where 
examples of our density profiles are plotted. Our grid of profiles, 
which are given names reflecting v x (e.g. "exp9" matches W7 at 
v x = 9000 km s -1 ), covers well the physical parameter space in 
which solutions for SN 2009dc will be found (see below). 

Table lA2l gives the factors a (for a given to) and slopes b of 
our empirical density profiles (cf. Equation IA2t . The value of b 
only depends on the rise time (assumed to be 21 d, 22.5 d, 25 d, 
27.5 d, 30 d, 35 d or 40 d). The parameter a depends on the as- 
sumed rise time as well, but furthermore on the choice of v x (9000, 
10000, 12000 or 140 00 km s" 1 ), at which the profile intersects W7 
jNomotoetaT]ll984h . We have chosen the reference time to to be 
31.5 d. 



8 For a hypothetical feature forming in a thin zone directly above the pho- 
tosphere, the curvature effect which makes the difference from line-of-sight 
velocities to radial velocity is most extreme: absorption is possible at line- 
of-sight velocities from km s up to the radial velocity. The radial veloc- 
ity then differs by a factor of 1.33 from the average line-of-sight velocity 
(calculated integrating over the lines of sight, assuming no limb darkening). 

9 W7 is evaluated for the reference time to we use in Equation iAH . 



Table A2. Parameters a and b of our empirical, exponential density models 
"exp..." for different choices of v x (exp9: v x = 9000 kms -1 , explO: 
v x = 10000 km s , ...) and assumed rise times t r . The a values are 
given for a reference time of to = 3 1 .5 d. 



assumed t r 


a (exp9 model) a 


(exp 10) 


a (exp 12) 


a fexpl4) 


b (all models) 


H 


[10- 12 gcm 


" 3 1 


[io- 12 g 


cm- 3 ] 


[10- 9 skm- 1 ] 


21.0 


15.39 


28.33 


77.21 


324.8 


8.25 


22.5 


8.93 


15.48 


37.38 


139.3 


7.64 


25.0 


4.21 


6.71 


13.70 


43.2 


6.81 


27.5 


2.29 


3.42 


6.10 


16.8 


6.13 


30.0 


1.39 


1 .96 


3.14 


7.7 


5.58 


35.0 


0.65 


Q.83 


1.12 


2.3 


4.72 


40.0 


0.37 


0.44 


0.53 


1.0 


4.09 



APPENDIX B: PARAMETERS OF THE ABUNDANCE 
TOMOGRAPHY MODELS 

The code input parameters for our Abundance Tomography models 
(Sections[3]|4]and[5]( are listed in Table IB ll Apart from the abun- 
dances, the code takes as input the photospheric velocity v p h, the 
time from explosion t and the bolometric luminosity Lboi- In addi- 
tion to the input values, Table IB ll also gives the calculated temper- 
ature of the photospheric black body, T p h, for each model. 
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Table Bl. Parameters of the models. The abundance values, which are only given for significant elements, apply between the photospheric velocity v p ^ of the 
respective epoch and of the previous epoch. The outer shells ("epoch" designations "oo", "o") improve the fit to the —9.4 d spectrum (cf. main text). For 
these shells, tipj, is not a photospheric velocity, but the velocity of the lower boundary of the respective zone. 



Model epochs lg ( L J ^lowbnd ^ph Element abundances (mass fractions) 





[d] 




Ikmr 1 ] 


[K] 


X(C) 


X(O) 


X(Mg) 


X(Si) 


X(S) 


A'(Ca) 


X(Ti) 


Jt(Cr) 


Jt(Fe) a » 


X( 56 Ni) a| 


09dc-AWD3det 


00 




11700 




0.03 


0.96 


0.00 


0.00 


0.00 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


t r = 30.0 d 







9650 




0.02 


0.S3 


0.04 


0.07 


0.04 


0.0000 


0.0000 


0.0001 


0.0000 


0.0008 




-9 


9.90 


91)50 


14264 


0.02 


0.56 


o.os 


0.15 


0.10 


0.0020 


0.0009 


0.0090 


0.0000 


0.0800 




-4 


10.01 


8300 


13931 


0.02 


0.45 


0.06 


0.25 


0.10 


0.0060 


0.0010 


0.0100 


0.0000 


0.1100 




+5 


9.98 


7900 


11519 


0.01 


0.27 


0.02 


0.45 


0.10 


0.0100 


0.0005 


0.0100 


0.0000 


0.1200 




+ 13 


9.82 


7500 


9224 


0.01 


0.15 


0.01 


0.50 


0.06 


0.0060 


0.0002 


0.0080 


0.0000 


0.2500 




+23 


9.65 


7200 


7591 


0.01 


0.08 


0.00 


0.50 


0.04 


0.0030 


0.0000 


0.0040 


0.0250 


0.5500 




+36 


9.39 


4700 


7291 


0.00 


0.01 


0.00 


0.05 


0.0 1 


0.0015 


0.0000 


0.0020 


0.2300 


0.7000 


09dc-AWD3det 


00 




11700 




0.03 


0.96 


0.00 


0.00 


0.00 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


E(B - V) = 0.06 mag 







9300 




0.02 


0.S9 


0.02 


0.05 


0.02 


0.0000 


0.0000 


0.0001 


0.0000 


0.0006 


t r = 27.8 d 


-9 


9.72 


8700 


14178 


0.02 


0.76 


0.04 


0.09 


0.03 


0.0015 


0.0006 


0.0055 


0.0000 


0.0550 




-4 


9.81 


8000 


13308 


0.02 


0.60 


0.04 


0.18 


0.06 


0.0080 


0.0008 


0.0080 


0.0000 


0.0800 




+5 


9.79 


7350 


11227 


0.02 


0.27 


0.04 


0.40 


0.10 


0.0080 


0.0015 


0.0150 


0.0000 


0.1500 




+ 13 


9.66 


7000 


9104 


0.01 


Q.18 


0.02 


0.45 


0.07 


0.0040 


0.0005 


0.0150 


0.0000 


0.2500 




+23 


9.5 1 


6650 


7518 


0.01 


0.14 


0.01 


0.30 


0.04 


0.0020 


0.0001 


0.0060 


0.0000 


0.5000 




+36 


9.26 


4200 


7346 


0.00 


0.1 1 


0.01 


0.25 


0.03 


0.0005 


0.0001 


0.0050 


0.0500 


0.5500 


uvuc-nvyujuct 




















n nnnn 
u.uuuu 


n nnnn 


n nnnn 




n nnnn 


solar progen. metallicity 


a 




10300 




0.02 


0.81 


0.03 


O.OS 


0.06 


0.0001 


0.0001 


0.0003 


0.0012 


0.0000 


t r = 30 d 


-9 


9.88 


9700 


13511 


0.02 


0.42 


0.06 


0.25 


0.10 


0.0053 


0.0001 


0.0030 


0.0650 


0.0750 




-4 


9.99 


8750 


13544 


0.02 


0.2S 


0.05 


0.36 


0.12 


0.0090 


0.0003 


0.0080 


0.0500 


0.1000 




+5 


9.97 


8300 


11098 


0.02 


0.02 


0.02 


0.65 


0.10 


0.0143 


0.0003 


0.0100 


0.0200 


0.1500 




+ 13 


9.82 


8000 


8921 


0.01 


0.01 


0.01 


0.60 


0.04 


0.0120 


0.0002 


0.0080 


0.0010 


0.3000 




+23 


9.66 


7800 


7306 


0.01 


0.01 


0.00 


0.40 


0.02 


0.0060 


0.0001 


0.0060 


0.0500 


0.5000 




+36 


l ).39 


4700 


7331 


0.00 


0.01 


0.00 


0.05 


0.00 


0.0010 


0.0000 


0.0030 


0.2300 


0.7000 


09dc-exp 






1 1700 




04 


95 


00 


00 


00 


6 0606 

1 1 1 11 I| |l 1 


n ooon 

1 I.I II 11 II 1 


6 0606 

1 1 1 11 I| |l 1 


n noon 
i i.i ii ii ii i 


n noon 


t r = 29.0 d 







9700 




0.02 


0.74 


0.06 


0.12 


0.06 


0.0000 


0.0000 


0.0001 


0.0000 


0.0008 




-9 


9.60 


9100 


14319 


0.02 


0.56 


0.06 


0.20 


O.OS 


0.0036 


0.0008 


0.0080 


0.0000 


0.0750 




-4 


10.01 


8550 


13670 


0.02 


0.39 


0.06 


0.35 


O.OS 


0.0060 


0.0010 


0.0100 


0.0000 


0.0900 




+5 


9.97 


7800 


11910 


0.02 


0.25 


0.03 


0.50 


0.06 


0.0060 


0.0005 


0.0100 


0.0000 


0.1250 




+13 


9.82 


7100 


9969 


0.01 


0.15 


0.02 


0.40 


0.05 


0.0045 


0.0003 


0.0100 


0.0000 


0.3500 




+23 


9.62 


6500 


8206 


0.00 


0.16 


0.01 


0.25 


0.02 


0.0015 


0.0001 


0.0075 


0.0250 


0.5250 




+36 


9.35 


3500 


9205 


0.00 


0.04 


0.00 


0.15 


0.01 


0.0008 


0.0001 


0.0050 


0.1000 


0.7000 


09dc-int 


oo 




11600 




0.08 


0.90 


0.01 


0.0 1 


0.00 


0.0000 


0.0000 


0.0000 


0.0005 


0.0001 


t r = 26.0 d 


O 




9200 




0.03 


0.54 


0.10 


0.25 


O.OS 


0.0010 


0.0001 


0.0001 


0.0015 


0.0030 




-9 


9.60 


8600 


14284 


0.03 


0.16 


0.10 


0.35 


0.15 


0.0025 


0.0010 


0.0100 


0.0500 


0.1500 




-4 


9.79 


7650 


14475 


0.03 


0.12 


0.04 


0.40 


0.15 


0.0050 


0.0010 


0.0100 


0.0005 


0.2500 




+5 


9.76 


7000 


11736 


0.03 


0.12 


0.02 


0.53 


0.15 


0.0050 


0.0008 


0.0080 


0.0005 


0.3500 




+ 13 


9.60 


6300 


9502 


0.01 


0.09 


0.01 


0.17 


0.05 


0.0050 


0.0004 


0.0070 


0.0500 


0.6000 




+23 


9.44 


5400 


8164 


0.00 


0.07 


0.01 


0.15 


0.01 


0.0050 


0.0002 


0.0070 


0.1000 


0.6500 




+36 


9.21 


4300 


7164 


0.00 


0.05 


0.00 


0.10 


0.01 


0.0050 


0.0001 


0.0070 


0.1500 


0.7000 



a * The abundances of Fe, Co and Ni in our models are assumed to be the sum of 56 Ni and the elements produced in the decay chain ( 56 Co and 56 Fe) on the one hand, and 
non-radioactive, directly-synthesised Fe (including also progenitor Fe) on the other hand. Thus, they are conveniently given in terms of the Ni mass fraction at t — 
[X( 56 Ni)o], the abundance of stable Fe at t — [X (Fe)nL ana " the time from explosion onset t. 
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